Homogenization of Elliptic Systems 
With Neumann Boundary Conditions 



Carlos E. Kenig* 



Fanghua Lin ' 



Zhongwei Shen* 



Abstract 



For a family of second order elliptic systems with rapidly oscillating periodic co- 
efficients in a C l,a domain, we establish uniform W 1,p estimates, Lipschitz estimates, 
and nontangential maximal function estimates on solutions with Neumann boundary 
conditions. 

1 Introduction and statement of main results 

The main purpose of this work is to study uniform regularity estimates for a family of ellip- 
tic operators {C £ ,e > 0}, arising in the theory of homogenization, with rapidly oscillating 
periodic coefficients. We establish sharp W l,p estimates, Lipschitz estimates, and nontan- 
gential maximal function estimates, which are uniform in the parameter e, on solutions with 
Neumann boundary conditions. 
Specifically, we consider 




(1.1) 



where e > 0. We assume that the coefficient matrix A{y) = (a^ (y)) with 1 < i,j < d and 
1 < a, (3 < m is real and satisfies the ellipticity condition 



m 2 < a$(y)&$ < -l^ 2 for y G R d and £ = (£f) e R dm , (1.2) 



where /i > 0, the periodicity condition 



A(y + z) = A(y) 



for y G M. d and z e Z d , 



(1.3) 



and the smoothness condition 



\A(x)-A(y)\ <r\x-y 



A 



for some A G (0, 1) and r > 0. 



(1.4) 
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We will say A G A(/i, A, r) if A = A(y) satisfies conditions (O]) . fOj) and fOD . 

Let / G L 2 (f2) and (7 G W~ l l 2,2 (dVt). Consider the Neumann boundary value problem 



£ £ «)=div(/) 
du £ 

g-n-f 



in Q, 
on dQ, 



where 



dv e _ 



ni(x)a"f(-) 



* e dxj 



;i.5) 



1.6) 



denotes the conormal derivative associated with C £ and n = {n\, . . . , n^) is the outward unit 
normal to dQ. Assume that j n u £ = 0. It is known from the theory of homogenization that 
under the assumptions fll.2l) - fll.3p . u £ — > u weakly in W 1,2 (Q) as e — > 0, where C (u ) = 
div(f) in f2 and |^ = g — n ■ f on <9f2. Moreover, the homogenized operator £ is an elliptic 
operator with constant coefficients satisfying (jl.2p and depending only on the matrix A (see 
e.g. 0). 

In this paper we shall be interested in sharp regularity estimates of u £ , which are uniform 
in the parameter e, assuming that the data are in LP or Besov or Holder spaces. The following 
three theorems are the main results of the paper. Note that the symmetry condition A* = A, 
i.e., 

a ij(y) = a ji"(y) f° r 1 — hi — d an d 1 < /3 < m, (1-7) 
is also imposed in Theorems 11.21 and 11.31 

Theorem 1.1 (W l,p estimates). Suppose A G A(yU, A,r) and 1 < p < 00. Let Q be a 

bounded C l > a domain for some < a < 1. Let g = (g p ) G B-V™(d£i), f = (/?) G L p (tt) 
and F = (F p ) G L q (Q), where q = ^ for p > and q > 1 for 1 < p < Then, if 

F and g satisfy the compatibility condition J n F^+ < g 13 , 1 >= for 1 < /3 < m, the weak 
solutions to 

( C £ (u £ ) =div(/) + F mil, 

dUs f AO 

g — n • j on oil, 



dv F 



(1.8) 



G W 1,p (f2) 



satisfy the estimate 



iVit 



e||lp(«) 



<C{ 



Z/P(n) 



|^||i9(Q) + IIS'IIb-i/p.p^)} 



1.9) 



where C > depends only on d, m, p, q, fi, X, r and Q. 



Theorem 1.2 (Lipschitz estimates). Suppose that A G A(/i, A,r) and A* = A. Let Q be a 

bounded C 1,a domain, < 77 < a < 1 and q > d. Then, for any g G C^dfi) and F G L 9 (£7) 
wt/i J Q F + g = ; the weak solutions to 



du £ 



dv F 



in Q, 
on dQ, 



;i.io) 



|Vw e | G 



2 



satisfy the estimate 

1 1 Vu e || i,oo(n) < C{\\g\\ C v(dn) + H-^IU^n)}, (1.11) 
where C > depends only on d, m, 7], q, fi, X, r and Q. 

Theorem 1.3 (Nontangential maximal function estimates). Suppose that A G A(/i, A,r) 
and A = A* . Let Q be a bounded C 1,a domain and 1 < p < oo. Then, for any g G L p (dQ) 
with mean value zero, the weak solutions to 

C £ (u £ ) = in O, 

< ^ = g on an, (i.i2) 

ov £ 

k (Vu £ y g L p (dn), 

satisfy the estimate 

||(Vu e )*|| LP (an) + ||VM £ || i9(n) < C ||p||LP(9n), (1.13) 
where q = jpr and C > depends only on d, m, p, /i, X, r and f2. 

A few remarks on notation are in order. In Theorem 11.11 B~ l l p,p (dVL) is the dual of the 
Besov space B l l p > p (dQ) on d£l, where p' = and < g@, 1 > denotes the action of g 13 on 
the function 1. By a weak solution u to ( 11.81) . we mean that u G W /1,p (fi) and satisfies 

L< (?) g ■ d £^= L {-^ + M d - < >■ (ll4) 

for any (p = (f a ) G Cq(M ). In Theorem 1 1.3 1 we have used (Vm £ )* to denote the nontangential 
maximal function of Vu E . We point out that the Lipschitz estimate in Theorem [T72] is sharp. 
Even with C°° data, one cannot expect higher order uniform estimates of u e , as V« E is known 
to converge to Viio only weakly. As a result, the use of nontangential maximal functions in 
Theorem 11.31 to describe the sharp regularity of solutions with L p Neumann data appears to 
be natural and necessary. Also note that under the conditions (11.21) and (11.41) . the existence 
and uniqueness (modulo additive constants) of solutions to (II. 8p . (jl.lOp and (11.12p with 
sharp regularity estimates are more or less well known (see e.g. [TJ [2J, [31]). What is new here 
is that with the additional periodicity assumption ( II. 3p . the constants C in the regularity 
estimates (11.91) . (II. lip and (11. 13j) are independent of s. 

In the case of the Dirichlet boundary condition u £ = g on <9f2 with g G B l ' p ' ,p {dVt) or 
g G C 1,v (dQ), results analogous to Theorems 11.11 and 11.21 were established by Avellaneda 
and Lin in [31 [7] for C l,a domains (without the assumption A* = A). They also obtained 
the nontangential maximal function estimate ||(tt e )*IU p (dn) < CH^Hlp^^) for solutions of 
C £ (u e ) = in Q (the case m = 1 was given in |3j). As it was noted in [3], uniform regularity 
estimates, in addition to being of independent interest, have applications to homogenization 
of boundary control of distributed systems [251 ESI E]- Furthermore, they can be used 
to estimate convergence rates of u e — > uq as e — > 0. In particular, it was proved in [3] 
that \\u £ — Wo|U°°(n) = 0(e), if C £ (u £ ) = div(/) in Q, u £ = g on dfl, and f,g are in 
certain function spaces. Extending the Lipschitz estimate ( II. lip to solutions with Neumann 
boundary conditions has been a longstanding open problem. The main reason why it is more 
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difficult to deal with solutions with Neumann boundary conditions in Theorem 11.21 than 
solutions with Dirichlet boundary conditions in [21 [7] is that now the boundary conditions 
in (11. 101) are e-dependent, which causes new difficulties in the estimation of the appropriate 
boundary correctors. We have overcome this difficulty, in the presence of symmetry, thanks 
to the Rellich estimates obtained in [211 122] . Neumann boundary conditions are important in 
applications of homogenization (see e.g. [SI [HI [261 12T])- The uniform estimates we establish 
in this paper can be used to study convergence problems for solutions u £ , eigenf unctions 
and eigenvalues with Neumann boundary conditions. As an example, let w e (x) = u e (x) — 
Uq(x) — ex{^)'Vuo(x), where x denotes the matrix of correctors for C £ in IR d . It can be 

shown that w £ = We + We , where || Vwe \\lp(u) < C p e|| V 2 mo||lp(q) for any 1 < p < oo, 
and \Vwi 2 \x)\dist(x } d£l) < Ce\\ VmoIIl 00 ^) f° r any £ G fi. We will return to this in a 
forthcoming publication. 

Let N £ (x,y) denote the matrix of Neumann functions for C e in Q (see Section 5). As a 
consequence of our uniform Holder and Lipschitz estimates, we obtain the following bounds, 

\N £ (x,y)\< ( 



d-2 ■ 



\x - y\ 

\V x N £ (x,y)\ + \V y N £ (x,y)\ < j^^TI' ( L15 ) 

C 



\V x V y Ne(x,y)\ < 



\x-y 



for d > 3 (see Section 8). In view of the work of Avellaneda and Lin on homogenization of 
Poisson's kernel [6], we remark that the techniques we develop in this paper may also be used 
to establish asymptotics of N e (x,y). This line of research, together with the convergence 
results mentioned above, will be developed in a forthcoming paper. 

We should mention that the case p = 2 in Theorem ll.3l is contained in [22J. In fact, for the 
elliptic system C £ (u £ ) = in a bounded Lipschitz domain Q, the Neumann problem with the 
uniform estimate || {Vu £ )*\\LP(dn) < C\\ §^ ||i>(8n) and the Dirichlet problem with the estimate 
II (u £ )*\\LP(dn) < C||w e ||Lp(en), as well as the so-called regularity problem with the estimate 
||(Vw e )*||rp(an) < C\\ VtanUe\\Lp(dfi), were solved recently by Kenig and Shen in [22] for p close 
to 2 (see [19] for references on boundary value problems in Lipschitz domains for elliptic 
equations with constant coefficients). The results in [22J are proved under the assumption 
that A G A(/i, A, r) and A* = A, by the method of layer potentials. In the case of a single 
equation (m = 1), the LP solvabilities of Neumann, Dirichlet and regularity problems in 
Lipschitz domains with uniform nontangential maximal function estimates were established 
in [21] for the sharp ranges of p's (the result for Dirichlet problem in Lipschitz domains was 
obtained earlier by B. Dahlberg [11], using a different approach; see the appendix to [21] for 
Dahlberg's proof). The results in [21~1 122] rely on uniform Rellich estimates || ^r\\L 2 (dn) ~ 
|| Vtan^e || l 2 (9c) f° r solutions of C £ {u £ ) = in a Lipschitz domain Q. We point out that 
one of the key steps in the proof of Theorem 11.21 uses the Rellich estimate || Vu £ \\ L 2^ dn ^ < 
c \\^\\L2(dn) in a crucial way. 

We now describe the key ideas in the proofs of our main results. To show Theorem II. 1[ 
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we first establish the uniform boundary Holder estimate for local solutions, 

\\ue\\c°^(B(Q,p)nn) < Cp' 1 (-1 \u £ \ 2 dx\ , (1.16) 

\JB(Q,2p)nn / 

for any 7 G (0, 1), where £ £ {u £ ) = in B{Q,3p) n ft and = on B{Q,3p) D <9ft for 
some Q € <9f2 and < p < c. The proof of (11.161) uses a compactness method, which was 
developed by Lin and Avellaneda in [21 El E] for homogenization problems, with basic ideas 
originating from the regularity theory in the calculus of variations and minimal surfaces. As 
in the case of Dirichlet boundary condition, boundary correctors are not needed for Holder 
estimates with Neumann boundary condition. From f l 1 . 1 6 [) one may deduce the weak reverse 
Holder inequality, 

(-f \Vu £ \ p dx) <C p (-f \Vu £ \ 2 dx) (1.17) 

for any p > 2. By [15] this implies that || Vw £ ||£p(q) < C|U ||i>(fi) for p > 2, if C £ {u £ ) = div(/) 
in Q and |^ = — n ■ f on d£l. The rest of Theorem 11.11 follows by some duality arguments. 

The proof of Theorem 11.21 is much more difficult than that of Theorem 11.11 Assume 
that G dQ. After a simple rescaling, the heart of matter here is to establish the uniform 
boundary Lipschitz estimate for local solutions, 

II Vw e ||Loo( B ( 0jl ) n Q) < C{ II u £ \\L^(B(o,2)nn) + ||fl'||c i '(B(o,2)nan)}, (1.18) 

for some i] > 0, where C £ (u £ ) = in -8(0, 3) Pi SI and ^f- = g on -8(0, 3) fl dVl. This problem 
has been open for more than 20 years, ever since the same estimate was established in [3] 
for local solutions with the Dirichlet boundary condition u £ = in -8(0, 3) fl dQ. Our proof 
of (11.181) also uses the compactness method mentioned above. However, as in the case of 
the Dirichlet boundary condition, one needs to introduce suitable boundary correctors in 
order to fully take advantage of the fact that solutions of the homogenized system are in 
C 1,r, (B(0, 2) D fi). A major technical breakthrough of this paper is the introduction and 
estimates of such correctors $ e = ($°^), where for each 1 < j < d and 1 < /3 < m, 
$f • = ($^-, • • • , ) is the solution to the Neumann problem 



£ £ (<,) = in n, 



^<*W = 011 3a (119) 



*^(o) = 0. 

Here Pj 3 = Xj(0, • • ■ ,1,...,0) with 1 in the (3 th position and J-^ denotes the conormal 
derivative of w associated with the homogenized operator Cq. Note that by the boundary 
Holder estimate, $"j(x) — > Xj5 a/ 3 uniformly in Q as e — > 0. To carry out an elaborate 
compactness scheme in a similar fashion to that in [3], one needs to prove the uniform 
Lipschitz estimate for the solution of (I1.19p . 

||V$ e |Uco (n) < C. (1.20) 
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The proof of (jl.20p relies on two crucial observations. First, one can use Rellich estimates 
as well as boundary Holder estimates to show that 



where \x — z\ < cdist(x, dQ). Secondly, if w £ (x) = § £ (x) — xL — ex(x/e), then ^f- can be 
represented as a sum of tangential derivatives of with ||<7y ||z°o(an) < Ce. Since C £ {w £ ) = 

in Q, it follows from these observations as well as interior estimates that |Vu> £ (x)| < 
C£[dist(a;, 90)] _1 . This gives the estimate |V$ £ (x)| < C, if dist(x,<9fi) > e. The remaining 
case dist(x, dQ) < e follows by a blow-up argument. See Section 7 for details. We note that 
the symmetry condition A* = A is only needed for using the Rellich estimates. 

With the Lipschitz estimate in Theorem 1 1.2 1 at our disposal, Theorem 1 1.31 for p > 2 follows 
from the case p = 2 (established in [22] for Lipschitz domains), by a real variable method 
originating in [9] and further developed in [28], [29J, [30]. The case 1 < p < 2 is handled by 
establishing L l estimate for solutions with boundary data in the Hardy space iJ 1 (9fi) and 
then interpolating it with L 2 estimates, as in the case of Laplacian [12] (see Section 9). In 
view of the Lipschitz estimates in [3] for local solutions with Dirichlet boundary condition 
and the L 2 estimates in [22], a similar approach also solves the LP regularity problem with 
the estimate || (Vu £ )*\\LP(dn) < C\\ VtanU £ \\LP{an) in a C 1,a domain Q for all 1 < p < oo (see 
Section 10). We further note that the same approach works equally well for the exterior 
domain Q_ = M d \ Q and gives the solvabilities of the LP Neumann and regularity problems 
in Q_. Consequently, as in the case of the Laplacian on a Lipschitz domain [321 E2], one may 
use the IP estimates in Q and f2_ and the method of layer potentials to show that solutions 
to the LP Neumann and regularity problems in C 1,Q domains may be represented by single 
layer potentials with density functions that are uniformly bounded in L p . Similarly, the 
solutions to the LP Dirichlet problem may be represented by double layer potentials with 
uniformly LP bounded density functions (see Section 11). 

The summation convention will be used throughout the paper. Finally we remark that 
we shall make little effort to distinguish vector- valued functions or function spaces from their 
real- valued counterparts. This should be clear from the context. 

2 Homogenization and weak convergence 

Let C £ = — div(A(x / e) V) with matrix A(y) satisfying fll.2|) - fll.3p . For each 1 < j < d and 

1 < /3 < m, let Xj — {Xji • • •■> Xj 1 ^) be the solution of the following cell problem: 



where Pj — Pj(y) — %(0, . . . , 1, . . . , 0) with 1 in the (3 th position. The matrix x = x(v) — 
{Xj{y)) with 1 < j < d and 1 < a, (5 < m is called the matrix of correctors for {C £ }. 




(1.21) 




(2.1) 
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With the summation convention the first equation in (12. ip may be written as 

d 



d_ 

dyi 



a0 

hj 



in 



(2.2) 



Let A 



where 1 <i,j < d, 1 < a, < m and 



-a/3 



[0,1]° 



a/3 i «7 



d_ 



xf 



dy. 



(2.3) 



Then £ = — div(Av) is the so-called homogenized operator associated with {C £ } (see [8]). 
We need the following homogenization result. 

Lemma 2.1. Let Q be a bounded Lipschitz domain in IR d and 

div [A k (x/e k ) Vu k ] =fe W~^ 2 (n) in SI, 

where e k — > and the matrix A k (y) satisfies ( f l.Hfy - fTlfy . Suppose that u k — > uq strongly in 
L 2 (fl), Vu k — > V«o weakly in L 2 (Q) and A k (x / e^Vuk converges weakly in L 2 (Q). Also 
assume that the constant matrix A k , defined by \2.3\) (with A replaced by A k ), converges to 
A . Then 

A k (x/e k ) Vu k — > A°Vu Q weakly in L 2 (fi) 
and div(A Vuo) — f in Q. 

Proof. If A k is independent of k, this is a classical result in the theory of homogenization 
(see e.g. [8J or [10]). The general case may be proved by the same energy method. We give 
a proof here for the sake of completeness. 



Let A k 



a/3 



l ij,k)' A k 



ijtk j and A = ). Suppose that 

dul 



weakly in L (f2). 



(2.4) 



Clearly, div(P) = / in fi, where P = (pf). For 1 < j,£ < d, 1 < /3 < m and k = 1, 2, 
write 

dul d 



dxe dx 
dx, a > h dx, 



£ kXTk( X / £ k) + X j$a/3 



' a ^' k dx~- { £kX Tk( X / £ k) + X 3 5 <*p} 



(2.5) 



where \* k — {Xjh) denotes the matrix of correctors for {C k )*, the adjoint operator of 
C k = — div(A k (x/e)'V). By taking the weak limits on the both sides of (12. 5p and using 
a compensated compactness argument (see e.g. Lemma 5.1 in [ID]), we obtain 



p,{x) ■ 



[0,1]° 



d_ 

dyi 



x7 k (y) + <M<*/3 }> dy 



OXp fc^oo 



[o,i] d 



07 

l il,k 



_d_ 

dyi 



xThiv) + <%<W \ d y- 



Since 



a 7,k(y)w-{x* j a k{y)}dy= / a^ k (y)— {xTliv)} d V 



'[0,l] d u Vi k " v </[0,l 

(see e.g. [HI p. 122]), it follows that 



_ o . lim ~/3 7 

dXi k^oo j£ ' k 

In view of (12. Ah this finishes the proof. □ 
Let if) : R^ 1 R be a C 1 - 00 function such that 

^(0) = |VV(0)| = and ||V^||c-o(R--i) < M , (2.6) 

where a G (0, 1) and M > will be fixed throughout the paper. For r > 0, let 

D{r) = D(r,ip) = {(x',x d ) G R d : \x'\ < r and ip{x') < x d < ip(x') + r}, 

D(r) = D(r,ip) = {(x',x d ) eM. d : \x'\ < r and ip(x') - r < x d < ip(x') + r}, (2.7) 

A(r) = A(r, V>) = {(x',^')) e Kd : M < r }- 

Lemma 2.2. Let {ip k } be a sequence of C 1,a ° functions satisfying A2.6\) . Suppose that ip k — > 
ipo in C l (\x'\ < r) and {\\Vk\\L 2 (D(r,ip k ))\ is bounded. Then there exist a subsequence, which 
we still denote by {vk}, and v G L 2 (D(r,ip )) such that v k — > v weakly in L 2 (Q) for any 
nccD(ri ). ' 

Proof. Let w k (x',x d ) = v k (x',x d + ip k (x')), defined on 

D(r, 0) = {(x',x d ) : \x'\ < r and < x d < r}. 

Since {w k } is bounded in L 2 (D(r, 0)), there exists a subsequence, which we still denote by 
{w k }, such that w k —> w weakly in L 2 (D(r, 0)). Let v (x',x d ) = w (x',x d — if) (x')). It is 
not hard to verify that v k — > vq weakly in L 2 (Q) if Q CC D(r, tpo). □ 

The following theorem plays an important role in our compactness argument for the 
Neumann problem. Note that (12.81) is the weak formulation of div(A k (x / £fc)V 'u k ) = in 
D(r,Vfc) and Sf = 9k on A(r,ip k ). 

Theorem 2.3. Let {A k (y)} be a sequence of matrices satisfying U.2\) - in~3\) and {ip k } a 
sequence of C 1,aa functions satisfying $2.6}) . Suppose that 

/ A k (x/e k )Wu k ■ \7(fdx = / g k -ipda (2.8) 

JD{r,i> k ) ^A(r,V> fc ) 
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for any ip G Cq(D(t, iftk)), where e k — > and 

\Wk\\w 1 > 2 (D{r,^ k )) + |bfc|U 2 (A(r,Vfc)) < C. (2.9) 

Then there exist subsequences of {ijj k }, {uk} and {gk}, which we still denote by the same 
notation, and a function ip satisfying \2.1$ , go G L 2 (A(r, ipo)), u$ G W 1,2 (D(r, ipo)), a 
constant matrix A such that 

ip k ip in C l (\x'\ < r), 

g k (x',tp k (x')) g (x',tp (x')) weakly in L 2 (\x'\ < r), (2.10) 
u k (x', x d - ^fe(^O) ~> u o(x', x d - i> (x')) strongly in L 2 (D(r, 0)), 



and 



/ A Vu ■ Vy? dx = / go-^pda (2.11) 



/or any y? G C*o(D(r, ^o))- Moreover, the matrix A , as the limit of a subsequence of {A k }, 
satisfies the condition M.2\) . 

Proof. We first note that H 2 . 101) follows directly from (12.91) by passing to subsequences. To 
prove (12. lip , we fix ip G Cl(D{r, ^o))- Clearly, if k is sufficiently large, ip G Cg(D(r, ipk))- It 
is also easy to check that 

/ g k -ifda ^ g -<pda. 

JA(r,ipk) ^A(r,-0 o ) 

By passing to a subsequence we may assume that A k — > A . Thus it suffices to show that 
/ A k (x/e k )Vu k ■ Vtpdx / A°Vu ■ Vy? dx. (2.12) 

In view of Lemma [2.21 we may assume that {u k }, Vu k , and A k (x / 'e^Vuk converge weakly 
in L 2 (Q) for any Q CC D(r,ip ). As a result, {u&J also converges strongly in L 2 (Q). 
Now, given any 5 > 0, we may choose a Lipschitz domain Q such that C D(r,ipo), 

\[ A°Vu - Vipdx\ < 5/3 (2.13) 

JD(r,i> )\n 

and 

A k (x/e k )Vu k ■ V(fdx\ < 5/3 (2.14) 



(2.15) 



>D(r,i> k )\n 

for sufficiently large. Thus ( 12.121) would follow if we can show that 

/ A k (x/e k )Vu k ■ Vipdx ->• / A°Vu ■ Vy? dx. 
Jn Jn 



This, however, is a direct consequence of Lemma [2.11 since div(A k (x / e k )Vu k ) = in Q by 
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We end this section with the uniform interior gradient estimate, established in [3] by 
Avellaneda and Lin, for solutions of C £ (u £ ) = 0. For a ball B = B(x,r) in IR d , we let 
pB = B(x, pr). We will use f- E f to denote t4- J e f, the average of / over E. 

Theorem 2.4. Let A e A(p, A,r). Suppose that C £ (u £ ) = in IB. Then 

sup|Vu e |<cQf |Vw £ | 2 da^ , (2.16) 
where C depends only on d, m, p ; X, r. 



3 Boundary Holder estimates 

The goal of this section is to establish uniform boundary Holder estimates for C £ under 
Neumann boundary condition. Throughout this section we assume that A e A(p, A, r). 

Theorem 3.1. Let Q be a bounded C 1,a ° domain. Let p > and 7 G (0, 1). Suppose that 
C £ (u £ ) = in B(Q, r) D and |^ = g on B(Q, r) D <9f2 /or some Q € <9f2 and < r < r . 
TTien 

sup |w e |<C<(-/ |w e | P dx) + p||fl'||l,»(B(Q I r)nan) > , (3-1) 

B(Q,r/2)nCl [ \JB(Q,r)nn J J 

and /or x, y G -B(Q, r/2) fl 

|u e (a;) - u e (y)| < C (— — — \ <(-f \u £ \ p dx] + p\\g\\L°°(B(Q, r )ndn) \ , (3.2) 



r J I \J B(Q,r)C\Vl J 

where ro > depends only on fl and C > on d, m, p ; X, r, p, 7 and fl. 

Let D(p,ip) and A(p,if)) be defined by (12. 7p . By a change of the coordinate system it 
will suffice to establish the following. 

Theorem 3.2. Let 7 e (0, 1). Suppose that C £ (u £ ) = in D(p) and |^ = g on A(p) /or 
some p > 0. T/ien /or any 1,1/6 D(p/2), 

I |\7 f / f \ 1/2 



\u £ (x) -u e (y)\ < C y — - — J U-j-^^\u e \J + p\\g\\L°°(A(p)) } , (3.3) 

where D(p) = D(p,ip), A(p) = D(p,ip), and C > depends only on d, m, p, A ; r, 7 and 
(a , M ) *n (EJ). 

The proof of Theorem 13.21 uses the compactness method developed in [3j El E] for homog- 
enization problems. We begin with the well known Cacciopoli's inequality, 

/ \Vu £ \ 2 dx< C - I \u £ \ 2 dx + Cp\\g\\ 2 LHA{p)) , (3.4) 

JD(sp) { l ~ s ) P JD{tp) 

du £ 
dv E 



where < s < t < 1, C £ {u £ ) = in D(p) and |^ = g on A(p). The periodicity of A is not 
needed here. 

For a function u defined on S, we will use (u)s (and f§) to denote its average over S. 
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Lemma 3.3. Fix (3 G (0, 1). There exist Eq > and 9 G (0, 1), depending only on d, m, \x, 
\, t, (3 and (a ,M ), such that 



\u 



e-mD(e)\ 2 <0 W , (3.5) 



D(fi) 



whenever e < e , C £ (u £ ) = in D(l), ^f- = g on A(l), 

IM|l°°(A(i)) < 1 and 4 \u £ - (u^)d(i)\ 2 < 1- 

Jd(i) 

Proof. Let Cq = — div(A°V), where A is a constant matrix satisfying (II. 2ft . Let 0' = 
(1 + /3)/2. By boundary Holder estimates for solutions of elliptic systems with constant 
coefficients, 

4 \w- {w) D{r )\ 2 < Cor 2 ?' for < r < -, (3.6) 

JD{r) 4 

whenever £ (w) = in D(l/2), g = g on A(l/2), 



||#|U°°(A(i/2)) < 1 and / \w\ < \D(1)\, (3.7) 

JD(X/2) 

where Co depends only on d, m, (3, \i and (ckq, A^o). 

Next we choose # G (0, 1/4) so small that 2Cq6 2/3 < 6 2/B . We shall show by contradiction 
that for this 9, there exists £o > 0, depending only on d, m, /i, A, r, /3 and (ao, M ), such 
that (13. 5p holds if < e < 6q and u e satisfies the conditions in Lemma [3.31 

To this end let's suppose that there exist sequences {sk}, {Ak}, {u £k }, {g k } and {ip k } 
such that e k — > 0, A k G A(/i, A, r), ^ satisfies (12. 6p . 

4\K) = mD k (l), 

du £k (3.8) 
- — = g k on A fc (l), 

||</fc|U«(A fc (l)) < 1, f K fc - K>D fc (i)| 2 < 1 (3-9) 

and 

/ K-(^) Dfc(e) | 2 >^, (3.10) 

where = — div(A k (x/e k )V) , D k {r) = D(r,ip k ) and A fc (r) = D(r,ip k ). By subtracting a 
constant we may assume that (u^)D k (i) = 0. Thus it follows from (I3.9P and the Cacciopoli's 
inequality (13 ,4p that the norm of u £k in W 1,2 (D k (l/2)) is uniformly bounded. In view of 
Theorem 12. 3| by passing to subsequences, we may assume that 

^ k -> ip in C\\x'\ < 1), 

g k (x',ip k (x')) -> ^o(x',^o(^)) weakly in L 2 (|x'| < 1), (3.11) 
u £k (x',x d -ip k (x')) -+uq{x',x a -iI)q{x')) strongly in L 2 (D(1/2,Q)), 
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and 



( div(A°Vw ) = in£>(l/2,Vo), 
(70 onA(l/2^o) 



i 



(3.12) 



where A is a constant matrix satisfying (11.21) . 
Using (13.111) one may verify that 



\D k (r)\ \D (r)\, \\g \\L°°{A(i,i, )) < h (u £k ) Dk(r) -> (u ) Do(r) 



and 



(3.13) 



D k (r) 



A>(1) 



\u £k ~ (u £k ) Dk(r) \ 2 ->• / |u - (no)flb(r)f 

for any r G (0, 1], where D (r) = D(r,ip ). It follows that 

kl 2 < i, 

) 

7 l M o - (uo)i 

J Do (6) 

In view of and fl3TTi|) we obtain Q 2fi < C 6 213 ' . This contradicts 2C 6 213 ' <Q 2fi . □ 

Lemma 3.4. Fix (3 G (0, 1). Let e , 6 be the constants given by Lemma \3lA Suppose that 
C £ {u £ ) = in D(l,vp) and §^ = g on A(l,ip). Then, if e < 6^% for some k>l, 



(3.14) 



)D ( 



> 



J \u £ -{u;) D{ek)i)) \ 2 <d 2 ^J 2 

J D{e k ,-4>) 



where 



J = max 



D(1,V0 



(3.15) 



K - (u £ ) D{m \ 2 ^j 



1/2 



, IMIl°°(A(i,v>)) 



Proof. The lemma is proved by induction on k. Note that the case k — 1 is given by Lemma 
13.31 Assume now that the lemma holds for some k > 1. Let e < ^^o- We apply Lemma [3.31 
to w(x) = u{6 k x) in D(l,tp k ), where ip k {x) = 6~ k ip{6 k x). Since C £ / k{w) = in D(l,ipf.)-i 
this gives 



7 l U £ ~ ( M e)D(e fe + 1 , 

J D(0''+ 1 

= 7 k-^W^fc)! 1 

JD(6,il> k ) 



max 



= ^max(-/ k-(^W^)| 2 , ^NIL) 
l</_D(6> fc ,v>) J 



where ||<?||oo = ||5 , ||l oo (A(i,V')) an d the last step follows by the induction assumption. Here we 
also have used the fact that || VipkWcoiRd- 1 ) < II V^Hc^o^-i) < M . □ 
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Proof of Theorem 13.21 By rescaling we may assume that p = 1. We may also assume 
that e < e , since the case e > e follows directly from the classical regularity theory. We 
may further assume that 



\uJ 2 < 1. 



IMU°°(A(i)) < 1 and / 

Jd{1) 

Under these assumptions we will show that 

\u £ - (T^) D (r)\ 2 < Cr 2(S (3.16) 



'D(r) 

for any r e (0, 1/4). The desired estimate (13. 3p with p = 2 follows from the interior estimate 
(I2.16P and (I3.16p . using Campanato's characterization of Holder spaces (see e.g. [IB"]). 

To prove (13.161) we first consider the case r > (e/eo). Choose k > so that 9 k+1 < r < 6 k . 
Then e < Eqt < Eo6 k . It follows from Lemma [3.41 that 



\U £ - (M £ ) D(r) | 2 < C 4 \U £ - {u e ) D{ek) \ 2 

D(r) J D(9 k ) 

< CQ 2kfi < Cr 2 ?. 

Next suppose that r < (e/e ). Let w(x) = u £ (ex). Then C\{w) = in D(eq , ip £ ), where 
ip e (x') = e~ l ij)(ex'). By the classical regularity we obtain 



\u £ - {u e ) D (r,i>)\ = 4 \w - (w) D (L^ s )\ 



\w - (w) D ,± M \ 2 , £ 2 \\g\ 



c [ 7) max \4 K- (u £ ) D (J-,i>)\ 2 , e 2 \\g\ 



00 



<C[-) (-1 =Ce^r 2 ? 



2/3 

where the last inequality follows from the previous case r = (e/eo). This finishes the proof 
of (I3.16P and thus of Theorem 13.21 □ 

We are now in a position to give the proof of Theorem 13.11 
Proof of Theorem 13.11 By rescaling we may assume that r = 1. The case p = 2 follows 
directly from Theorem 13.21 To handle the case < p < 2, we note that by a simple covering 
argument, estimate (13. ip for p = 2 gives 

sup \u £ \ <C< 1 (-f lu.l 2 ] + \\g\\ L °°(B(Q,i)ndn) \ , (3.17) 
B(Q,s)nn [ l J - s ) \JB(Q,t)nn J J 

where (1/4) < s < t < 1. By a convexity argument (see e.g. [T4J p. 173]), estimate (I3.17P 
implies that for any p > 0, 

[4 K\ 2 ) <cA(4 \u e \ p +\\ g \\ L (3.18) 

\J B{Q,i/2)nn J [ \JB(Q,i)nn J J 
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The case < p < 2 now follows from estimate (I3.18P and the case p = 2. □ 

4 Proof of Theorem 11.11 

Under conditions fll.2j) and ( II. 4p . weak solutions to (I1.8P exist and are unique, up to an 
additive constant, provided that the data satisfy the necessary condition J Q F l3 + < g@,l >= 
for 1 < (3 < m. In this section we will show that the weak solutions satisfy the uniform 
W 1,p estimate in Theorem ll.il 

Our starting point is the following theorem established by J. Geng in [15], using a real 
variable method originating in [9] and further developed in [281 1211 130] . 

Theorem 4.1. Let p > 2 and Q be a bounded Lipschitz domain. Let £ = — div(A(:r)V) be 
an elliptic operator with coefficients satisfying (£Q|). Suppose that 

>i/p r r ■) 1/2 

\vu\n <Co\4 \vu\ 2 } , (4.i) 



Bnn ) KJ2Bnn 



whenever ue W 1 ' 2 {3Bnfl,) ! C{u) = m35n0, and gf = on 3BndQ. Here B = B(Q,r) 
is a ball with the property that < r < r and either Q 6 dQ or B(Q,3r) C il. Then, for 
any f e L P (Q), the unique (up to constants) W 1 ' 2 solution to 

C(u) = div(/) in Q, 

— — = — n ■ f on ml, 
av 

satisfies the estimate 

||V«||iP (n) < C p \\f\\ LP(n) , (4.3) 
where C p depends only on d, m, p, fi, ro, Q and the constant Cq in (f^._/|). 



Now, given A e A(/x, A, r) and p > 2. Let Q be a C 1 ' 00 domain. Suppose that C £ {u £ ) = 
in 3B n Q and |^ = on 3B n <9fi. If 35 C fi, the weak reverse Holder inequality (jCT]) for 
m £ follows from the interior estimate ( I2.16p . Suppose that Q G dfl and = B(Q,r). We 
may use the interior estimate and boundary Holder estimate (13.21) to obtain 



\Vu £ (x)\ < C5(xY l (J \u £ (y) - u £ (x)\ 2 dy] 

\Jb(x,c5(x)) J 



6(x) 



<C, ttt [4 Vu £ \ 2 dy 



1/2 

2 



B(Q,2r)nn 



1/2 

(4.4) 



for any 7 G (0, 1) and x G B(Q,r) fl fi, where 5(x) = dist(x, dfl). Choose 7 G (0, 1) so that 
P7 < 1. It is easy to see that (14. 4 p implies 



f / \Vu e A ^ < C p f-f \Vu £ \ 2 ) ' 
\JBnn J \J2Bnn / 

In view of Theorem 14. II we have proved Theorem II. II for the case p > 2, g = and F — 0. 
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Lemma 4.2. Suppose A G A(/x, A,r). Lei / G L p (fi) ; where Q be a bounded C l,a ° domain 
and 1 < p < oo. Let it G be a weak solution to C £ (u £ ) = div(/) in Q and |^ = — n-f 

on dQ. Then \\Vu £ \\ L p(q) < C p \\f\\ L p(n)- 

Proof. The case p > 2 was proved above. Suppose that 1 < p < 2. Let g G C^°(Q) and v £ 
be a weak solution of £e(v e ) = div(g) and = on <9f2, where £* denotes the adjoint of 
C £ . Since A* G A(A, /i,r) and > 2, we have ||Vu £ || iP '(Q\ < C'Hfi'll^'^)- Also, note that 



SI C^i J \£/ ctej * dXi 

where / = (/") and g = (<7f). The estimate || Vu e ||i>(n) < C||/||LP(n) now follows from (14. 5 \ 
by duality. □ 

Lemma 4.3. Suppose that A G A(A,/z,r). Lei g = (g Q ) G B~ 1 ^ p ' p (dft), where Q is a 
bounded C 1,a ° domain, 1 < p < oo and < g Q , 1 >= 0. Let u G W 1,P (Q) be a weak solution to 
C £ {u £ ) = in SI and ^ = g on dtt. Then \\Vu £ \\ LP{n) < C p \\g\\ B -Vp,p{da)- 

Proof. Let / G C£°(tt) and v £ be a weak solution to C*(v e ) = div(f) in f2 and |^f- = on 
Since A* G A(A,/x, r), by Lemma l4~2l we have ||Vv e ||_ z y^ < C H/Hjym) 1 
Note that 

Let £" be the average of t> £ over fi. Then 

| < fiS^e > | = | < g,v £ - E > I < \\g\\ B -i/P'P(dn)\\ v e - E\\ B i/ p , p '^ 
— C \\g\\B- 1 /p-p(an)\\ v e — E\\ w i, p '^ 

< C lbllB- 1 /p.P(9n)||Vf e ||2 /P '(n) 

< C||5 , llB- 1 /p i p(aQ)ll/llLp'(n)5 



(4.7) 



where we have used a trace theorem for the second inequality and Poincare inequality for the 
third. The estimate ||Vw e ||LP(n) < C \\g\\B- 1 /P'P(an) follows from (l4.6p -( ]4~Tj) by duality. □ 

Let 1 < q < d and \ = \ ~\- ^ n ^ ne P ro °f °f ^ ne next lemma, we will need the following 
Sobolev inequality 

\ 1 /p f r \ V? 

«dx < C / |Vw| 9 cfe , (4.8) 



where u G VF 1,9 (f2) and j m u = 0. 

Lemma 4.4. Suppose that A G A(/i, A,r). Lei F G L 9 (fi) ; where 1 < q < d and Q is a 
bounded C 1 ' 010 domain. Let u G W l,p (VL) be a weak solution to C £ (u £ ) = F inVl and ^f- = —b 
on dQ, where \ = \~\ and b = ]aJT\ In F - Tfien \\^ u e\\Lp(n) < C \\F\\ Lq ^. 
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Proof. Let / G C£°(fi) and v £ be a weak solution to (C £ )*(v £ ) = div(/) in Vl and |^ = on 
<9f2. By Lemma [4.21 we have || Vi^U^'m) < C 11/11 lp'(q)- Note that 

8U * f?dx= I afi-YS-^dx 



q dx-i J q ^ V s / dxj dx 



F-v £ dx- b-v £ da (4.9) 



F(w £ - £) da;, 

where E is the average of v £ over <9f2. It follows from (14. 9 ft and Sobolev inequality ( 14.81) that 



r du a 

j n ^i:-f? dx \<\\ F \\L< m \\v £ -E\ 



Li' (Si) 



< C||i r IU9(n)||Vu e || i y( n ) 



< C||F||L9(n)||j || LP ' (n) . 

By duality this gives ||Vw e ||iJ.(n) < C ||F||i 9 (n). □ 

Proof of Theorem ll.il Let v £ be a weak solution to C £ (v £ ) = div(/) in f2 and |^ = —n-f 
on <9fi. Let u> £ be a weak solution to C £ (w £ ) = F in f2 and = —6 on dCl, where 
b = J n F. Finally, let h £ = u £ — v £ — w £ . Then C £ (h £ ) = in and |^ = g + 6 on <9fi. 
It follows from Lemmas I4.2[ 14.31 and 14.41 that 

||Vu e ||Li>(fi) < ||Vu £ ||iF(n) + \\Vw £ \\ LP (n) + \\Vh E \\ L p(si) 

< C {||/||ip(n) + + \\g\\B- 1 /p-p(dn)} , 

where q = for p > -^j, and q > 1 for 1 < p < -^j. This completes the proof. □ 

5 A matrix of Neumann functions 

Let T £ (x,y) = (r^ e (x,y)) denote the matrix of fundamental solutions of C £ in M, d , with 
pole at y. Under the assumption A e A(/i, A, r), one may use the interior estimate (I2.16P to 
show that for d > 3, 

\T £ (x,y)\<C\x-y\ 2 - d (5.1) 

and 

\V x T £ {x,y)\ + \V y T £ (x,y)\ < C\x-y\ l ~ d , (5.2) 

where C depends only on d, m, /i, A and r (see e.g. [T7]; the size estimate 05. ip also 
follows from [13]). Let V E (x,y) = (V^ £ (x,y)) mxm , where for each y e Q, Vf(x,y) = 

i V Ae( x > 2/)> • • • > V^f 2/)) solves 

C £ (Vf(-,y))=0 infi, 

k -im Jdsi 
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where T^(x, y) = (r A y), . . . , T™ e (x, y)) and e 13 = (0, . . . , 1, . . . , 0) with 1 in the f3 
position. We now define 

N £ (x,y) = (N^ £ (x,y)) mxm = T £ (x,y) - V e {x,y), 

for x,y ett. Note that, if Ng(x, y) = V%(x, y) - Vf(x, y), 



th 



(5.4) 



C e {Nfi(;y)} = ef'8 v (x) in ft, 
£-{Nl(',y)} = -eP\Xl\- 1 ondQ 



(5.5) 



N?(x,y) da{x) = 0, 



mi 



where 6 y (x) denotes the Dirac delta function with pole at y. We will call N e (x, y) the matrix 
of Neumann functions for C F in ft. 



(5.6) 



Lemma 5.1. For any x,y G ft, we have 

where A* denotes the adjoint of A. 
Proof. Note that 

T Ae( x ' V) = r Av(l/> x )' for an y y G ft. (5.7) 
Using the Green's representation formula for C £ on ft, (15.31) and (15. 5p one may show that 

1 



Vg{x,y)+T%{x,y) 



a/3 



ajf ( - ) n» Jz,x 



z\ d 



e) dz. 



This gives V^ £ (x,y) = V^ £ (y,x) and hence (15.61) . 



□ 



Theorem 5.2. Let Q be a bounded C 1,a ° domain and A G A(p,X,r). Lei XojJ/oj^o 

G ft 6e 

sitc/j i/iai |xq — z \ < (1/4) [rro — yo\- Then for any 7 G (0, 1), 



(/ \V y {N £ (x ,y) - N E (z ,y)}fdy) 1/2 < Op 1 -" ()*B-*L) 

U B( yo , P /4)nn J \ P J 

where p = \xq — yo\ and C depends only on p, X, r, 7 and ft. 

Proof. Let / G C$°(B(y , p/2) n ft) and J n f = 0. Let 



(5i 



u e (z) = / N £ (x,y)f(y)dy. 
Jo. 
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Then C £ (u £ ) = f in Q and = on dtt. Since C £ (u £ ) = in B(x Q ,p/2) n ft, it follows 



from the boundary Holder estimate (13. 3p and interior estimates that 



\u s (x ) - u £ (z ) \ < C 



Fo — z o 



\Vu F 



P J VJ B(x , P /2)nn 

Let E be the average of u £ over S(y , p/2) fl f2. Note that by ( II. 2p . 



1/2 



(5.9) 



/ \Vu £ \ 2 dx < 


/ f -u £ dx 




L 


In 


Jq 







f ■ (u £ - E) dx 



B(yo,p/2)ntt 

< ||/IU 2 (n)|K - E\\ L 2 ( B (y ,p/2)nn) 

< Cp\\f\\ L 2(ty\\Vu E \\ L 2( B (y 0)P / 2 )nn), 



(5.10) 



where we have used the Cauchy and Poincare inequalities. Hence, || Vu £ \\ L 2^ < Cp\\f\\ L 2my 
This, together with (15. 9p . gives 



\u £ (x ) - U £ (z )\ <Cp 2 a ' 



L 2 (Q)- 



By duality this implies that 



-| 1/2 

W(y)-C X0)Z0 pdy^ <Cp 

B(y ,p/2)nU ) 



2,1' _ 2-1 / FO — 2 o| 



P 



(5.11) 



where H^(|/) = N £ (x ,y) — N £ (z Q ,y) and is the average of W over B(y ,p/2) fl fi. In 

view of (ESD we have (C e )*(W*) = in B(y ,p/2) n ft and ^{W 7 *} = on dfi, where ^ 
denote the conormal derivative associated with (£ e )*. The estimate ( 15 .8p now follows from 
( 15. lip by Cacciopoli's inequality ( 13 .4p . □ 



Lemma 5.3. Let V £ (x,y) be defined by ( 15. 3\) . Suppose d >3. Then for any x,y E Q, 

\V £ (x,y)\<C[6(x)]^[6(y)] 2 -^, (5.12) 

where 5(x) = dist(x, dfl). 

Proof. We begin by fixing y e Q and 1 < f5 < m. Let u £ (x) = V £ (x,y). In view of (15. 3p we 
have 

du £ du £ 
||Vu £ || L 2( n ) < — \\ w -i/2,2 {dn) < C|| — || iP (an), 



dv F 



where p = Note that by (E2]) 

„ du £ „ 



< C 
<C[5(y 



2-d 
2 



Thus we have proved that 



\Vu £ \\ L 2 (n) <C[5(y)\ 
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Now, by the interior estimates and the Sobolev inequality (14.81) . 
\u £ (x)\ < C \ [ \u £ {z)f dz\ ' 

<C[5(x)]^\ ( f \Vu £ \ 2 dx] +\n\&\-[ u £ da\ 



2-d 



<C[5(x)] » |[%)] » +\n\^\-f T e (z,y)da(z)\ 
<C[5(x)]^[5(y)]^, 



mi 



where 2* = j^. □ 

Theorem 5.4. Let Q be a bounded C 1,a domain in M. d , d > 3. Suppose that A e A(/i, A, r). 
Then 

\N £ (x,y)\<C\x-y\ 2 - d (5.13) 

and /or any 7 6 (0, 1), 

\N E (x,y)-N E (z,y)\< C " 



|JV 6 (y,a;)-JV e (y >iZ )|< r ^_ 

F — 2/1 



_ id-2+7' 



-2+7 ; 



where \x — z\ < (1/4) |x — y|. 

Proof. By Theorem 13.11 we only need to establish the size estimate (I5.13p . To this end we 
first note that by Lemma [5.31 

\N £ (x,y)\ < C {\x-y\ 2 - d +[5(x)] 2 - d +[5(y)] 2 - d }. (5.15) 

Next, let p — \x — y\. It follows from Theorem 13 . 1 1 and (I5.15P that 

i/p 

\N s (z,y)\*dz\ 

B(x,p/4)nn 
<C{\x-yt d + [5(y)] 2 - d }, 



{£ 



\N e {x,y)\<CU4 \N £ {z,y)\ p dz} + 



\dn\ J (5.I6) 



where we have chosen p so that p(d — 2) < 1. With estimate (I5.16P at our disposal, another 
application of Theorem 13.11 gives 



\X:U.U)\<C{U \N e (x,z)\>dz\ P -\ ,r- 1 

U B(y,p/4)nSl 



<C\x-y\ 2 - d . 

This finishes the proof. □ 
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Remark 5.5. If m = 1 and d > 3, the size estimate (1 5 . 1 3 [) and Holder estimate (I5.14p 
for some 7 > were established in [2U] for divergence form elliptic operators with bounded 
measurable coefficients in bounded star-like Lipschitz domains. 

Remark 5.6. Suppose that d > 3. The matrix of Neumann functions for the exterior 
domain = R d \ U may be constructed in a similar fashion. Indeed, let N £ (x, y) = 
T £ (x, y) — V~(x, y), where V~(x, y) is chosen so that for each y G 0_, 

' £ £ {N~(-,y)} = 6 y (x)I inQ, 

< JL{N-(-,y)} = ondn, (5.17) 

N~(x,y) = 0(\x - y\ 2 ~ d ) as |x| ->■ 00, 

where I is the m x m identity matrix. The estimates in Theorem 15.41 continue to hold for 
N-(x,y). 

Remark 5.7. If d = 2, the matrix of Neumann functions may be defined as follows. Choose 
B(0,R) such that Q C B(0,R/2). Let G e (x,y) be the Green's function for C £ in B(0,R). 
Define N £ (x, y) = G £ (x, y) — V E (x, y), where V e (x, y) is the solution to (I5.3p . but with V £ (x, y) 
replaced by G £ (x, y). Theorem [5]2] continues to hold for d = 2. One may modify the argument 
in the proof of Lemma 15.31 to show that 

\V e (x,y)\ <C,[6(x)]-"[6(y)]-\ 

for any 7 > 0. In view of the proof of Theorem 15.41 and the estimate \G £ (x,y)\ < C{1 + 
I In \x — y\\\ in [3], this gives \N £ (x,y)\ < C 7 |x — y|~ 7 for any 7 > 0. 



6 Correctors for Neumann boundary conditions 

Let $ £ = where for each 1 < j < d and 1 < < m, = ($JJ, . . . , $ is a 

solution to the Neumann problem 

im^is?!) (6 - i) 

Here Pj 3 = Pj(x) — Xj(0, . . . ,1, . . . , 0) with 1 in the (3 th position. In the study of boundary 
estimates for Neumann boundary conditions, the function $ e (x) — x plays a similar role 
as sx{~) f° r interior estimates. The goal of this section is to prove the following uniform 
Lipschitz estimate of $ £ . 

Theorem 6.1. Let Q be a C 1,a ° domain. Suppose that A e A(/i, A, r) and A* = A. Then 

||W £ || L co (n) < C, (6.2) 
where C depends only on d, m, \i, \, r and Q. 
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Our proof of Theorem 16.11 uses the uniform L 2 Rellich estimate for Neumann problem: 

\Vu £ \ 2 da <C [ \^\ 2 da, (6.3) 



an Jan 



for solutions of C £ (u £ ) = in Q. We mention that (16. 3p as well as the uniform L 2 Rellich 
estimate for the regularity of Dirichlet problem: 



\Vu £ \ 2 da< C / \V Ua u £ \ da, (6.4) 
an Jan 

was established by Kenig and Shen in [22] under the assumption that Q is Lipschitz, A G 
A(/i, A, r) and A = A* (also see [21] for the case of the elliptic equation). The constant C in 
(I6.3p - fl6.4l) depends only on d, m, /x, A, r and the Lipschitz character of Q. 



Lemma 6.2. Let Q and C satisfy the same assumptions as in Theorem \6.1[ Suppose that 



C £ {u £ ) — in Q, |^ = g on d£l, and 



where gij G C x (dQ) and n = (m, • • • , w<j) denotes the unit outward normal to dQ. Then 

\Vu E (x)\ < ^rX)||<7til|L»(an), (6-5) 
/or any x G fi, where 8(x) = dist(x, dQ). 

Proof. By the interior estimate (I2.16P we only need to show that 

\ue(x) -u £ (z)\ < C^Hflfylliaofan), (6-6) 

where \x — z\ < cr and r = o~(x). Let N £ (x,y) denote the matrix of Neumann functions for 
C £ on Q. Note that 

u e (x)-u B (z)= / {N £ (x,y) - N E (z,y)}g(y)da(y) 
Jan 

= ~ S ^lj d ~ nj ^y~) i Ns ( x >y} ~ N ^ Z ^V)} ■ 9ij(y)do-(y), 

where we have used the fact that tangential derivative on dVt. Consequently 

it suffices to show that 

/ \V y {N £ (x,y)-N £ (z,y)}\d<j(y)<C, (6.7) 
Jan 
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if | a; — z\ < cr and r = 5(x). 

Let Q G <9fi so that \x — Q\ = dist(x, dQ). By translation and rotation we may assume 
that Q = and 

Q fl : | a;' | < 8cr and |x<z| < 8cr} 

= Ux'jXd) : \x'\ < 8cr and ip(x') < Xd < 8cr} 
where ip(0) = |V^(0)| = and c is sufficiently small. To establish (16. 7p we will show that 



/ |V„{iV e (a:,y)- N e {z,y)}\da(y) < C, 

J \y\<cr 

and there exists f3 > such that for cr < p < r , 

[ \V y {N £ (x,y) - N £ (z,y)}\dcr(y) <c(-)' 

J\y-P\<cp \PJ 



(6.8) 



(6.9) 



where P G dQ and \P\ = p. The estimate (16.71) follows from (16.81) and (16. 9p by a simple 
covering argument. 
To see (JSSD we let 

S'(t) = {(x',X(t) : | a;' | < t and ^(x') < < ^(x') + ct}. 

Note that by Cauchy inequality, for t G (cr, 2cr), 



\y\<cr 



\V y {N £ (x,y) - N £ (z,y)}\da(y) 



<Cr d - 1 [ \V y {N £ (x,y)-N £ (z,y)}\ 2 da(y) 

JdS(t) 



< Cr 



d-l 



dS(t) 



(6.10) 







du': 



{N £ (x,y)-N £ (z,y)} 



da(y), 



where we have used the Rellich estimate ( 16. 3 p for the last inequality. Since 

g^{N £ (x,y)-N £ (z,y)}=0 in dSl, 
we may integrate both sides of ( I6.10P in t over (cr, 2cr) to obtain 



y\<cr 



\V y {N £ (x,y) - N £ (z,y)}\da(y) 



<Cr d - 2 [ \V y {N £ (x,y)-N £ (z,y)}\ 2 dy. 

JS{2cr) 



(6.11) 



The desired estimate ( 16. 8 p now follows from estimate ( 15. 8p . 
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The proof of (16. 9p is similar to that of (16. 8p . Indeed, an analogous argument gives 
( / \V y {N £ (x,y)-N £ (z,y)}\ 2 da(y)\ 

VJ\y-P\<cp ) 

<Cp d - 2 [ \V y {N £ (x,y)-N £ (z,y)}\ 2 dy 

J\v-P\<2cp 



\y-P\<2cp 
27 

< c ' ' 



This completes the proof. □ 
Let ^ £ = where 1 < j < d, 1 < a, (3 < m and 

= - - exf (f ) • (6.12) 

Lemma 6.3. Suppose that fl and £ satisfy the same conditions as in Theorem \6 . 1[ Then 

Ce 

|V* e (x)| < —, — r for any igSI. (6.13) 
5{x) 

Proof. Fix 1 < i < d and 1 < 7 < m. Let w = {w\ . . . , w m ) = . . . , ^™ 7 ). Note that 

C £ {w) = in f2. In view of Lemma 16.21 it suffices to show that there exists G C 1 (dQ) 
such that 

dw \ s f d d 



dv £ 



(6.14) 



To this end we observe that by the definition of $"j in (I6.ip . 

/ 9W \ n a ( ' X\ 8 f _ a^~\ n a (X\ 8 



I ^— I = n,a°f J */»T I 



" 1 ij \e) dxj \ £/ J \e) dxj 



~cry a/3 
Wit ~ Wj 



3l 



; 7 



dxj V £ , 

where a*? are the homogenized coefficients defined by (12.31) . Let 



H^(V) = a7 - agiy) \ + ^(V) } ■ (6-15) 



It follows from the definition of that 



H?{y)dy = Q. 
[o,i] d 
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Thus we may solve the Poisson equation on [0, l] d with periodic boundary conditions, 



r ° 7 in R d , 



f AU? = El 
1 U°[ f (y) is periodic with respect to Z d . 



(6.16) 



Since A(y) and V%(y) are Holder continuous, V 2 U°P is Holder continuous. In particular, we 
have llVL^Iloo < C, where C depends only on /i, A and r. 
Now let 



dyk 



Then 



and hence 



cry 



rn(x) 



d 

dx k 



x 



P«7 
tr ilk 



We claim that 



Assume the claim is true. We may then write 







x 



on dVL. 



(6.17) 



(6.18) 



(6.19) 



Since ||ei^(a?/e)||oo < C £ > we obtain the desired f !6.14|) . 
Finally, to show ( I6.18p . we observe that 



d_ 

dyt 



H7(V) = o 



in 



which follows directly from (12.21) . In view of (16.161) this implies that ^p{U°f(y)} is harmonic 
in W 1 . Since it is also periodic, we may deduce that -^{U^iy)} is constant. As a result, 



d_ 

dyi 



Fuliv) 



d 2 



dy k dy, 

This completes the proof of Lemma 16. 2[ 

Proof of Theorem EZQ It follows from ( 1021) and lETI5j) that 



in 



|v*.(,)| < c , 



for any x E CI. 



□ 



(6.20) 



This implies that |V<E> e (x)| < C if 5(x) > c£. To estimate |V$ e (x)| for a; with S(x) < ce, we 
use a standard blow-up argument. 
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Fix j and f3. Let w(x) = e A £x )- Then C\{w) = and 

tt- = -T^-yex) = n i {ex)a i -. 
ov\ av £ J 

Since fl is a C 1,a ° domain, its normal n(x) is Holder continuous. Thus, by the classical 
regularity results for the Neumann problem with data in Holder spaces, 

||V$ e |U« (B( Q iB)n n) <C + c\\ [ \V$ £ \ p dx) ^ (6.21) 

l £ JB(Q,2e)nn J 

for any p > 0, where Q G dQ and C depends only on d, m, p, p,, A, r and fl. We remark 
that estimate (I6.2ip with p = 2 is well known and the case < p < 2 follows from the case 
p = 2 by a convexity argument. Finally, it follows from (I6.20p and (I6.2ip with p < 1 that 

||V (: I ) e ||L°°( J B(Q,£)nQ) < C- 

This finishes the proof of Theorem 16.11 □ 

Remark 6.4. Fix rj G C^(R d ' x ) so that r){x') = 1 for \x'\ < 2 and r](x') = for \x'\ > 3. 
For any function ip satisfying the condition ( 12. 6p . we may construct a bounded C l,a ° domain 
Q$ in M. d with the following property, 

D^(i) Gfl^ G {{x',x d ) : \x'\ < 8 and \x d \ < 8(M + 1)}, 
{(x',(^)(x')) : \x'\ <4} Cffi f 

Clearly, the domain can be constructed in such a way that fl^ \ {(x r , (iprj)(x r )) : Isc'l < 4} 
depends only on M . 

Let <& e (x) = $> £ (x,Q^, A) be the matrix of functions satisfying ( 16. ip with Q = and 
$ e (0) = 0. It follows from Theorem 16.11 that ||V$ e || £«>($}) < C, where C depends only on d, 
m, p,, A, t and (a ,M ). 



7 Boundary Lipschitz estimates 

In this section we establish the uniform boundary Lipschitz estimate under the assumption 
that A G A(/x, A, r) and A* = A. 

Theorem 7.1. Let Q be a bounded C l,a ° domain. Suppose that A G A(/i, A, r) and A* = A. 

Let C £ (u £ ) = in B(Q, p) fl Vl and ^ = g on B(Q, p) fl dfl for some Q G dfl and < p < c. 
Assume that g G C V (B(Q, p) fl dVL) for some r\ G (0,«o)- Then 

||Vn e || LO o(£(Q )P /2)nn) < C {p~ ||« e ||L°°(B(Q,p)nn) + ||5 , ||c;(B(Q,p)n9C) } , (7.1) 

where c = c(Q) > and C depends only on d, m, p, X, t, rj and fl. 

Let D(p) = D(p,iji) and A(p) = D(p,ip) be defined by (EZJ) with ijj G C^R^ 1 ), 
^(0) = |V^(0)| =0 and ||V^|| c a o(R d ' 1 ) < Mq. We will use HgHc ^^) to denote 

inf {M : \g(x) - g(y)\ < M\x - y\ p for all x,y G if}. 
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Lemma 7.2. Let < rj < «o and k = (1/4)77. = ^e{ x i^"4>iA) be defined as in 

Remark 6.4 There exist constants Eq > 0, 9 G (0, 1) and Cq > 0, depending only on d, m, 
ji, X, r, 7] and (a ,M ), such that 



\u £ - < $ e , B e > lU-p^) < 9 1+K , (7.2) 



for some B e = (b e ■) G R with the property that 



whenever 



|B e | < Co^lKIU.-^)) and < n(Q)A,B e >= n^affb^ = 0, 



£ < e , £ £ (m £ ) = m D(l), : — - = # on A(l), w e (0) = 0, 



and 

Ibllco.^ACi)) < 1, 9(0) = 0, ||u e || L oo (D{1)) < 1. (7.3) 

Proof. Let Co = — div(v4°V), where A = (fi£- ) is a constant m x m matrix satisfying ( II. 2p . 
By boundary Holder estimates for gradients of solutions to elliptic systems with constant 
coefficients in C l,a ° domains, 



\\W- < X, (Vw) D(r ) > ||L°°(D(r)) ^ ^ 

< Cir 1+2h - {||o||c(A(l/2)) + ||w||l°=(D(1/2))} , 

for any r G (0, 1/4), whenever C (w) = in D(l/2), g = g on A(l/2) and w(0) = 0. The 
constant Ci in (17.41) depends only on d, m, [A, rj and (ao, Mo). Observe that if 

o (0) =<n(0)A°,(V«;)(0) >= 0, 

then ||g||c"j(A(i/2)) < C||5 , ||c ."(A(i/2)) and 



< n(0)A u , (Vw) fl(r) > I = I < n(0)A u , (Vw) D(r) - (Vic)(0) > 

< Cr 2K {\\g\\c°>V(A(l/2)) + IMU°°(D(l/2))} • 



(7.5) 



Consequently, if we let B = (&o G R dm with 



l g^-.#M»)?(|) jH (7.6) 

where {h al3 ) mxm is the inverse matrix of (nj(0)nj(0)ay ) mXm , then 

< x, B > ||r«(iJ(o,r)) < ^ 2 r 1+2K , (7.7) 
for any r G (0, 1/4), provided that C {w) = in D(l/2), g = g on A(l/2), w(0) = 0, 

\\g\\c°>HA(i/2))<l> 0(0) =0 and |MU°°(£>(i/2)) < 1, (7-8) 
where C 2 depends only on d, m, /1, 77 and (a , M ). 
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Next we choose 9 G (0,1/4) so small that 2C2# K < 1. We shall show by contradiction 
that for this 9, there exists e > 0, depending only on d, m, //, A, r, r\ and (a ,M ), such 
that estimate ( I7.2p holds with 



if < e < 6q and u £ satisfies the conditions in Lemma [7.21 We recall that (af?) in (I7.9P is 

the homogenized matrix given by ( 12. 31) . It is easy to verify that rij(0)a"- fofj = 0. Also, by 
the divergence theorem, |B e | < Co0 -1 ||« £ ||i,oc,m(0)). 

To show ( I7.2p by contradiction, let's suppose that there exist sequences {e k }, {A k }, {u ek }, 
{g k } and ipk such that e k — > 0, A k G A(/i, A, r), satisfies ( 12. 6p . 



f <K) = in2? fc (l), 

^ = 9k on A fc (l), (7.10) 
u £k (0) = g k (0) = 0, 

bfc||co."(A fc (l)) < 1, ||« Efc |U°°(O fc (l)) < i, (7.ii) 



and 



K fc -<<,B>|| L ^ (e)) ># 1+K , (7.12) 

where D fc (r) = D(r,ip k ), A k (r) = A(r,ip k ), $^ = § £k {x,VL^ k , A k ) and is given by (EH]). 
By passing to subsequences we may assume that as k — > oo, 

A k ->• A , 

V>*->$> in^dx'l <4), (7.13) 
g k {x',ij; k (x')) ->■ ^(a 7 , V'oOO) in ^(k'l < !)■ 

Since ||« £ Jc^CD(i/2,v fc )) + ||*e fc ||c7(D(i/2,-0 fc )) < C by TheoremO, again by passing to subse- 
quences, we may also assume that 



u £k (x',x d - ip k {x')) -> M (x',x d - ^o(^')) uniformly on £>( 1/2,0), 
— ipk(x')) converges uniformly on D(l/2, 0), 

where R^ Ax) = Q k k (x) — x. Furthermore, in view of Theorem 12.31 we may assume that 
£oM = in D(l/2,Vb) and ^ = g on A(l/2,z/>o), where £ = -div(A°V). 

Note that by Lemma |6\3| R^ Ax', Xd — ipki.%')) must converge to a constant. Since (0) = 
0, we deduce that R k k (x',Xd — ip k (x')) converges uniformly to on .D(l/2, 0). Thus, in view 
of (I7~n|) - (l7~np . we may conclude that u (0) = #(0) = 0, 

IM|c°."(A(l/2,V>o)) < 1) || m o||l°°(D(1/2,Vo)) - 1 (7-15) 

and 

||«o- < x, B > |U»(£,(^)) > 9 1+K . (7.16) 

This, however, contradicts with ( I7.7p -( I7T8|) . □ 
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Remark 7.3. Let w =< $ e ,B e >= ^(z)^ where $ £ and B e are given by Lemma 17721 
Then C £ {w) = and J^- = rii{x)a^b^. In particular, we have w(0) = and J^r(O) = 0. 
Also, note that in Lemma |7.2[ one may choose any 9 e (0,#i), where IC-fl* = 1. These 
observations are important to the proof of the next lemma. 



Lemma 7.4. Lei « ; eo? # be the constants given by Lemma 7.2. Suppose that C £ (u £ ) = in 
D(l,ip), ^f- = g on A(l,tp) and u £ (0) = g(0) = 0. Assume that e < 9 1 ~ 1 e§ for some £ > 1. 
Then there exist B{ e R dm for j = 0, 1, . . . , i - 1, such that 



< n(0)A,B £ >= 0, |B e | <CJ 



and 



where 



e-i 

\u £ -J2 6Kj < U l> B i > < (7-17) 

j=0 



n £ (x)=^(r^,^,A), 

J = max {||5'||c '"(A(i,v))7 \\u £ \\l°°{D(i,^)) } 

andipj(x') = 9- j ?p(9 j x'). 

Proof. The lemma is proved by an induction argument on £. The case I — 1 follows by 
applying Lemma 17.21 to u £ /J. Suppose now that Lemma 17.41 holds for some £ > 1. Let 
£ < 9 e eo. Consider the function 



w(x) 



e~ £ <j u £ {9 e x) - 0Kj < n i(#H b{ > 

j=0 



on D(l,i[)e). Note that Cs r {w) = in D(l,ipt), w(0) = and by the induction assumption, 

(9* 

\\w\\ L o°(d(im) < 9 £k J. (7.18) 



Let 



Then 



dw 



/i(a;) = £(# £ x) - ^ < n{9 e x)A, B £ >, (7.19) 
i=i 

where n denotes the unit outward normal to A(l,^>). It follows that h(0) = 0. Since 
s6 < Eq, we may then apply the estimate for the case £ = 1 to obtain 

\\W~ < $£(x,tyfc,A),B£ > 

„ 1+K r ii 7 ii ii n i v ' 

< max {||/i||c .i(A(i,^)) ; IM|i/>°(d(i,^)) / , 

where Be e IR dm satisfies the conditions < n(0) A B e >= and 

} • (7-21) 
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It follows that 

e-i 

< 9 i+1+K max { ||/i||c°'')(A(i,^)) ) IIHU 00 ^!,^))} • 
To estimate the right hand side of (17.221) . we observe that 

C°.«»(A(1,V><)) 



j=0 



1 C||^llc ."(A(l,V>)) 

since 77 = 4k. Since < 77 < ao 5 by making an initial dilation of x, if necessary, we may 
assume that |M|c ' , '(A(i,i/>)) ^ s sma U so that 

g<t | 1 + ^IN^(A(i,v)) |< L (723) 

This implies that 

||^||c°.»j(A(i,v><)) ^ ^ (7.24) 
This, together with flTTTSl) and f l7T22l) . gives 

|| tte _ J20 Kj < ni > 11^(^+1^)) < fl^W+'O J, (7.25) 
3=0 

where we have chosen = 8~ £k ~B « . Finally, in view of (I7.2ip . ("IT. 18[) and (I7.24p . we have 

(9* 

|B^| < CJ. This completes the induction argument. □ 
Lemma 7.5. Suppose that C £ (u £ ) = in D(l) and |^ = g on A(l). TTien 

/ |Vtt e | 2 dar< Cp d {IMioop^ + lbH^i))}, (7-26) 

/or any < p < (1/2), where C depends only on p, A, r, 77 and (M ,a ). 

Proof. By subtracting a constant we may assume that u £ (0) = 0. We may also assume that 
g(Q) = 0. To see this, consider 

where (bP) G M m solves the linear system nj(0)nj(0)a°^6 /3 = g a (0). Then C £ {v £ ) = in D(l), 
v £ (0) = and 

29 



Thus |^(0) = 0. Since ||$ e ||i,° o (D(i)) + ||V$ e ||i,oo(£)(i)) < C, the desired estimate for u £ follows 
from the corresponding estimate for v £ . 

Under the assumption that u e (0) = g(0) = 0, we will show that 

\\u E || l°°(d( p )) < Cp {||w £ ||l<x>(d(i)) + IM|c"j(A(i))} , (7.27) 

for any < p < (1/2). Estimate (I7.26P follows from ( I7.27P by Cacciopoli's inequality (13. 4p . 

Let k, So, 9 be the constants given by Lemma 17.21 Let < e < 9e$ (the case e > 9e$ 
follows from the classical regularity estimates). Suppose that 

9 i+1 < — < 9 i for some i > 1. 

Let p G (0,1/2). We first consider the case — < p < 9. Then 6 e+1 < p < p l for some 
i = 1, . . . , %. It follows that 



(7.28) 



|pe|U°°(.D(p)) < ll M e||L°°(D(0O) 

e-i £-1 

< \\ Ue -j2& Kj <niM>\\L°°iD { et)) +J2 dKj \ B i\\\ n i\\^(D 

3=0 j=0 

e-i 

< j + C J J2 8 Kj Pill r« [Dm > 

3=0 

where J = max {Hs'llcMp^)), ||w £ ||.l°°(.d(i))} and we have used Lemma EH Recall that 
Ili(x) = 9 j $^(9~ j x, Qy j: A). By Remark EH we have IP £ (0) = and || VIF £ || L oo (D(1)) < C. 
Hence, 

l|n^||2,oop(0*)) < co . 

This, together with ( 17.281) . gives ||w £ ||l°° (£>(/>)) < CpJ for any — < p < | (the case 9 < p < 
(1/2) is trivial). 

To treat the case < p < — , we use a blow-up argument. Let w(x) = e~ 1 u E (ex). Then 
Ci(w) = in D(2sq 1 ,iI) £ ) and §^(x) = g(ex) on A(2£q 1 ,^ £ ), where ip e (x') = e~ 1 '^{ex'). By 
the classical regularity estimate, 

i „ '-(iii I "9w 



It follows that 

||Vu e || L oo( D (j^)) < C |£ _1 ||n|| L oo (D (|£)) + ||#||c»j(a(i))} < CJ, 

where we have used the estimate (17.271) with p = ^ for the last inequality. Finally, since 
w e (0) = 0, for < p < — , we obtain 



\u F 



\l°°[d[p)) < Cp|| Vm £ ||l°°(d(^)) < CW- 



This completes the proof of (17.271) . □ 
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Proof of Theorem 17.11 By rescaling we may assume that p = 1. By a change of the 
coordinate system, we may deduce from Lemma 17.51 that if P G dQ, \P — Q\ < ~ and 

< r < i, 

/ \Vu £ \ 2 dx < Cr d | ||' u e||Loo(_B(Q I i)nf2) + llfi , llc , i(B(Q,i)n9n)} > 

JB(P,r)nn 

where C depends only on d, m, /x, A, r, 77 and fl This, together with the interior estimate 
f )2.16p . implies that 

ll^ W e||l,°°(.B(Q,§)nfi) < C {l|We|U°=(B(Q,l)nn) + Hfl 1 II c(s(Q,i)n9n) } • 

The proof of Theorem 17.11 is now complete. □ 



8 Proof of Theorem 1.2 



Under the condition A G A(A, /x, r), we have proved in Section 5 that 

\N £ (x,y)\<- C ^- j - 2 if d>3. (8.1) 

With the additional assumption A* = A, we may use Theorem 17. II to show that for d > 3, 

{X -J l (8.2) 

iv.v^^y)! < i - d . 

\x — y\ a 



1-7 



If d = 2, one obtains |iV e (x, , y)\ < C 7 |x— y|~ 7 and |Va;vV e (x, y)\ + \V y N £ (x, y)\ < C^\x—y 
for any 7 > (this is not sharp, but sufficient for the proof of Theorem II .2p . Now, given 
F G L q (il) for some q > d, let 

v e {x) = / N £ (x,y)F(y)dy. 
Jn 

Then C £ (v £ ) = F in f2 and |^ = — J^i* 1 on <9f2. Furthermore, it follows from pointwise 
estimates on \V x N £ (x,y)\ that || Vv £ \\L°°(n) < C||-^||M(n)- Thus, by subtracting v £ from u e , 
we may assume that F = in Theorem 11.21 In this case we may deduce from Theorems 17.11 
and EO that for Q G dfi, 

||VM £ ||Loo( B (Q iP / 2 )nn) < C < I 4 |Vm £ | 2 ) + \\g\\cv(A(Q, P )) } , (8.3) 

[ \JB(Q,p)nn / J 

where C depends only on d, m, /x, A, r, 77 and Q. Since ||Vix e ||i2(n) < C||fi , ||L 2 (9n), the 
estimate || Vu £ \\l^(u) < C\\g\\cv(dn) follows from (18. 3 p and the interior estimate ( 12. 16ft by a 
covering argument. 
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9 Proof of Theorem 11.3 



As we mentioned in Section 1, the case p = 2 is proved in [22] (for Lipschitz domains). To 
handle the case p > 2, we need the following weak reverse Holder inequality. 

Lemma 9.1. Let fl be a bounded C 1,a ° domain. Suppose that A G A(A,/i,r) and A* = A. 
Then, for Q G dfl and < r < r , 

sup (Vu £ )*<c\-f \(Vu £ )*\ 2 do\ , (9.1) 



where u £ G W 1,2 (B(Q, 3r) D fl) is a weak solution to C £ (u £ ) = in B(Q, 3r) fl wi/i either 
|h* = o or u £ = on B(Q, 3r) n dfl. 

Proof. Recall that the nontangential maximal function of (Vu £ )* is defined by 
(Vtt e )*(P) = sup {\Vu £ (x)\ : x G Q and \x - P\ < C dist(x, dfl)}, 
for P G <9f2, where Co = C(fl) > 1 is sufficiently large. Note that 

(Vu £ )*(P) = m&x{Mr,l{Vu £ ),M r ,2{Vu £ )}, 

where 

M r ,i(Vu £ )(P) = sup {|Vm £ (x)| : x efl, \x — P\ < c r and \x — P\ < C dist(x, dfl)}, 
M r ^iyu £ )(P) = sup {|Vu e (a;)| : x G fl, |x — P| > c r and |x - P| < C dist(x, <9fi)}, 

and Co = c(fl) > is sufficiently small. Using interior estimate (I2.16p . it is easy to see that 
su PB(Q,r)ndn Mr^fVue) is bounded by the right hand side of (19.10 . To estimate JA r> i(Vu e ), 
we observe that 



sup A4 r ,i(V« £ ) < sup |Vw 

B(Q,r)ndCl B(Q,3r/2)nfl 



; I 



<C {4 \Vu £ \ z dx 

' B(Q,2r)nU 



1/2 

(9.2) 

1/2 

<C {4 \(Vu £ )*\ 2 da 

B(Q,2r)ndQ 

We point out that the second inequality in (19.21) follows from the boundary Lipschitz estimate. 
For Neumann condition |^ = on B(Q,3r) fl fl, the estimate was given by Theorem 17.11 
while the case of Dirichlet condition follows from Theorem 2 in [3] p. 805]. □ 

Lemma 9.2. Suppose that A G A(A,/i, r) and A* = A. Let p > 2 and fl be a bounded 
Lipschitz domain. Assume that 

\Vp / r \l/2 

\(Vu £ y\ p da) <C[4 \(Vu £ )*\ 2 da) , (9.3) 

B(Q,r)ndU J \JB(Q,2r)ndn J 

whenever u £ G W 1 ' 2 (B(Q, 3r) fl fl) is a weak solution to C £ {u £ ) = in B(Q,3r) fl fl and 
|^ = on B(Q,3r) fl dfl for some Q G dfl and < r < r$. Then the weak solutions to 
C £ (u £ ) = in fl and |^ = 36 L p (dfl) satisfy the estimate \\(Vu £ )*\\LP(dn) < C \\g\\LP(dn)- 

32 



Proof. This follows by a real variable argument originating in [9] and further developed in 
[2"8"| I2"9~| I3D] . In [21] the argument was used to prove that for any given p > 2 and Lipschitz 
domain Q, the solvability of the Neumann problem for Laplace's equation Am = in fi 
with L p boundary data is equivalent to a weak reverse Holder inequality, similar to (19. 3p . 
With the solvability of the L 2 Neumann problem for C £ (u £ ) = [22] and interior estimate 
(I2.16p . the proof of the sufficiency of the weak reverse Holder inequality in [24"| pp. 1819-1821] 
extends directly to the present case. We omit the details. □ 

It follows from Lemmas 19.11 and 19.21 that Theorem 11.31 holds for p > 2. To handle the 
case 1 < p < 2, as in the case of Laplace's equation [12J, one considers the solutions of the 
L? Neumann problem with atomic data |^ = a, where J dQ a = 0, supp(a) C B(Q,r) fl dQ 
for some Q G dQ and < r < r , and HaHi^a^) < r 1_d . One needs to show that 



The case 1 < p < 2 follows from (19. 4p by interpolation. 

To prove (19.41) . one first uses the Holder inequality and the L? estimate ||(Vw £ )*||t2( a Q\ < 

C \\a\\L 2 (dn) < Cr 1 ^ to see that 




(9.4) 




(9.5) 



Next, to estimate (Vw)* on dVt \ B(Q, Cr), we show that 




(9.6) 



for some 7 > 0, where p 



P -Q\>Cr. Note that 




Q)}a{y) da(y) 



(9.7) 



for some b G R m . It follows that 




Q)}\da(y) 



(9.8) 



Hence, if z G Q and cp < \z — P\ < Cq5(z) for some P G B(P , cp) fl dQ, 
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where 5(z) = dist(z, dQ) and we have used the interior estimate, Minkowski's inequality and 
Theorem 15.21 This implies that 



L 



x T 



-M 2 , p (Vm £ ) da<C[-\ . (9.9) 

B(P ,cp)n9n VP/ 

Finally, to estimate A^i,p(Vw e ), we note that the L 2 nontangential maximal function 
estimate, together with an integration argument, gives 



[ \M hp (Vu £ )\ 2 da<- [ \Vu £ \ 2 dx, 

J B(p ,cp)ndn P JB(p ,2cp)nn 

(see p2] for the case of Laplace's equation). It follows by Holder inequality that 

1/2 



(9.10) 



M hp {Vu £ ) da < Cp d - 1 (-1 \Vu £ \ 2 dx 

\J B(P ,2cp)nQ 



B(p ,cp)ndn 

< C 



(9.11) 



where the last inequality follows from (19 .8p and Theorem 15.21 In view of (I9.9P and (19. lip , 
we have proved (19.51) . The desired estimate 



/ (Vu £ )*da<C 

Jdn\B(Q,Cr) 



>dn\B(Q,Cr) 

follows from (19. 5 p by a simple covering argument. This completes the proof of (19. 4p and 
hence of Theorem 11.31 □ 

Remark 9.3. The estimate || Vu £ \\LQ(n) < C||^||i P (gn) with q = -J^j- in Theorem 11.31 follows 

from Theorem ll.il using the fact that L p (dVl) C B~~i ,q (dVt). The estimate also follows from 
the observation that ||u>||L<7(n) < C\\ (w)*\\lp(oq) for any w in a Lipschitz domain Q. To see 
this, we note that 



Jan \ x ~ vl 



By a duality argument, it then suffices to show that the operator 



*i(/)(*)= / ^rw^ dy 

Jn \ x ~ y\ 



is bounded from L q ' (Vl) to L p ' (dVl). This may be proved by using fractional and singular 
integral estimates (see e.g. [29, p. 712]). 

Remark 9.4. Suppose that d > 3. For g G L p (dQ), consider the L p Neumann problem in 
the exterior domain f2_ = IR d \ Q, 

C £ (u £ ) = in 

g on dn, (9.13) 



\2-d 



k (Vu £ )* G U\dQ) and u £ (x) = 0(\x\ 2 - d ) as |x| ->• oo. 
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It follows from [22] that if p = 2 and Q is a bounded Lipschitz domain with connected 
boundary, the unique solution to ( I9.13P satisfies the estimate ||(Vw £ )*||i2j-an) < C \\g\\L 2 (m) 
(if dQ is not connected, the data g needs to satisfy some compatibility conditions). An 
careful inspection of Theorem 11.31 shows that the L 2 results extend to LP for 1 < p < oo, if 
Q is a bounded C 1,a domain. 



10 L p Regularity problem 

In this section we outline the proof of the following. 

Theorem 10.1. Suppose that A G A(/i, A, r) and A* = A. Let Q be a bounded C 1,a domain 
with connected boundary and 1 < p < oo. Then, for any f G W 1,p {dVl), the unique solution 
to C £ {u £ ) = in Q, u £ = f on dfl and (Vu £ )* G L p (dQ) satisfies the estimate 

||(Vu e )1lp(«i) < C llVton/ll^sn), (10-1) 

where C depends only on d, m, p, \i, X, r and Q. 

The case p = 2 was proved in |22j for Lipschitz domains. The case p > 2 follows from 
Lemma 19.11 and the following analog of Lemma 19. 2[ 

Lemma 10.2. Suppose that A G A(A,/x, r) and A* = A. Let p > 2 and Q be a bounded 
Lipschitz domain with connected boundary. Assume that 

(7 \(Vu £ )Tda) ^ <C f-f |(Vn £ )*| 2 ^ ' , (10.2) 

\JB(Q,r)ndn / \JB(Q,2r)ndn / 

whenever u £ G W 1,2 {B{Q 1 3r)nfi) is a weak solution to C £ (u £ ) = in B(Q, 3r)nQ and u £ = 
on B(Q, 3r) fl dQ for some Q G dQ and < r < ro- T/ien f/ie weaA; solution to C £ (u £ ) = 
m Q and u £ = f G H /1,p (9fi) satisfies the estimate \\(Vu e )*\\LP(dn) < C ||Vt on /||jy»(an)j where 
C depends only on d, m, p, fi, X, r, r , Cq and Q. 

The proof of Lemma 110.21 is similar to that of Lemma 19. 2[ We refer the reader to [23] 
where a similar statement was proved for elliptic equations with constant coefficients. 

To handle the case 1 < p < 2, we follow the approach for Laplace's equation in Lipschitz 
domains [12] and consider L 2 solutions with Dirichlet data u £ = a, where supp(a) C B(Q, r)D 
dQ for some Q G dQ and < r < r , and || V fan a||ioo( 9n ) < r l ~ d . By interpolation it suffices 
to show estimate (19. 4p . Note that \a\ < Cr 2 ~ d . Using the estimates on Green's functions in 
[3], one has 

Cr 

|Vu E (a;)|< 1 ii\x-Q\>Cr. (10.3) 

\x — Q\ a 

Estimate f)9.4p follows easily from the L 2 estimate ||(V« e )*||2,a(am < C||Vton a IU 2 (8S2) an d 

Remark 10.3. One may also consider the LP regularity problem for the exterior domain: 
given / G W 1,p (dQ), find a solution u £ to C £ {u £ ) = in Q_ such that u £ = f on dfl, 
(Vm £ )* G L p (dtt) and u £ {x) = 0{\x\ 2 ~ d ) as |ar| -> oo. It follows from [22] that if is 
a bounded Lipschitz domain in R d , d > 3, then the unique solution to the L 2 regularity 
problem in Q_ satisfies the estimate ||(Vu e )*||ia(9n) ^ C \\^tanf\\w 1 ' 2 (dfi)- An inspection of 
Theorem 110.11 shows that the L 2 result extends to LP for 1 < p < oo, if Q is a C 1,a domain. 
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11 Representation by layer potentials 

For / G L p (dQ), the single layer potential u £ = S £ (f) and double layer potential w £ = T> £ (f) 
for the operator C £ in Q are defined by 



/ \ p (11-1) 

<(*) = / ( ^{rSvCv.*)} J f'WMv), 



e 



where r^ ^x, ?/) and TA*, e (x,y) = (r^ >e (y,x))* are the fundamental solutions for C £ and 
(£ £ )* respectively. Both S e (f) and V e (f) are solutions of £ e (tt) = in R d \ dQ. Under the 
assumptions that A G A(/i, A, r) and Q is a bounded Lipschitz domain, it was proved in [22] 
that for 1 < p < oo, 



(V5 e (/))l^ (n) + ||(P e (/)) < C7 f 



Lp(dn), 



where C p depends only on d, m, fi, A, r, p and the Lipschitz character of Q. Furthermore, 

(Vu e )±(P) exists for a.e. P G dQ, (§^) ± = (±|/+/C A>£ )(/) and K) ± = (^I+)C% jE )(f), 

where JC* A * £ is the adjoint of /C^* j£ . Here (w)± denotes the nontangential limits on dQ of u, 
taken from Q and f2_ respectively. 

Let L^(dQ,R m ) denote the space of functions in L p (dQ,R m ) with mean value zero. 

Theorem 11.1. Let Q be a bounded C l,a domain in R d , d > 3 with connected boundary. 
Suppose that A G A(/i, A, r) and A* = A. Then, for 1 < p < oo, 

l -I + TC A ,e ■■ L p (dQ,R m ) -+ L p (dQ,R m ), 

--I + )C% je : L p (dQ, R m ) ->■ L p (<9ft, M m ), ( 1L2 ) 

S £ : L p (<9fi,R m ) -»■ H /1,p ((9fi,lR m ), 

are invertible and the operator norms of their inverses are bounded by a constant independent 
ofe. 

Proof. The case p = 2 was proved in [22] for Lipschitz domains. If Q is C 1,a , the results for 
p 7^ 2 follow from the solvabilities of the L p Neumann and regularity problems with uniform 
estimates in Q and f2_ (see Theorem II. 3[ Theorem 110 .1| Remarks 19.41 and I10.3j) . □ 

As a corollary, solutions to the LP Dirichlet, Neumann and regularity problems for 
C £ {u £ ) = may be represented by layer potentials with uniformly LP bounded density 
functions. This shows that the classical method of integral equations applies to the elliptic 
system C £ {u £ ) = 0. 

Theorem 11.2. Let 1 < p < oo. Under the same assumptions on A and Q as in Theorem 
the following holds. 

(i) For g G L p (dQ), the solution to the LP Dirichlet problem in Q with u £ = g on dQ is given 
by u £ = V £ (h £ ) with ||/t e ||z,P(an) < C p \\g\\ LP ( d n). 
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(ii) For g E L p {dVt), the solution to the L p Neumann problem in f2 with |^ = g on dVl is 
given by u £ = S £ (h £ ) with \\h £ \\ LP ( dn ) < C p \\g\\ L p idn) . 

(Hi) For g e W 1,p (dQ), the solution to the L p regularity problem in Q with u £ = g on dfl is 
given by u £ = S £ (h £ ) with \\h £ \\ LP(dn) < C p \\g\\ LP(dQ) . 
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